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We have determined spectroscopic ages of elliptical galaxies in the Virgo clus- 
ter using spectra of very high signal-to-noise ratio (S/N A^^ >100). We observed 
8 galaxies with the Subaru Telescope and have combined this sample with 6 galax- 
Ph| ies previously observed with the WHT. To determine their ages we have used a 

O I new method based on the H70- age indicator, which is virtually independent of 

^ [ the effects of metallicity. Apart from ages we have estimated abundances of var- 

' ious elements. In this paper we present the observations, the data reduction and 
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1. Introduction 

An important issue in current-day astronomy is whether the stellar populations in el- 
liptical galaxies are uniformly old or a complex mixture of various populations of different 
ages. One of the classical formation scenarios, the monolithic collapse scenario (Eggen, 
Lynden-Bell & Sandage 1962) suggests that all eUiptical galaxies consist of uniformly old 
stars : i.e., they formed at high redshift {zf >2) and evolved passively. The fundamental 
plane (Djorgovski & Davis 1987; Dressier et al. 1987; Bender, Burstein & Faber 1993) and 
the Mg2 — a relation (Bender et al. 1993; Ziegler & Bender 1997; CoUess et al. 1999) show 
indications towards this possibility. On the other hand, the hierarchical clustering scenario 
predicts that elliptical galaxies were formed as a result of merging of galaxies accompanied 
by intensive starbursts (e.g. Toomre & Toomre 1972; Kauffman & Chariot 1998). Results 
based on spectroscopic line strength analysis by several authors (e.g. Gonzalez 1993; Trager 
et al. 2000b) seem to support this scenario, since a significant scatter in ages is found. Re- 
cent results from the SLOAN survey (Kauffmann et al. 2003) seem to indicate that large 
ellipticals are old and small ellipticals are younger. Although measuring accurate ages and 
metallicities would lead to a breakthrough in our understanding of the formation of elliptical 
galaxies, this has been considerably difficult due to the age-metallicity degeneracy (Worthey 
1994; Arimoto 1996). Photometric and spectroscopic properties such as colors, absorption 
line indices, change almost in the same way to differences in age and metallicity of the stellar 
populations, so that it is difficult to determine unambiguous ages. 

In understanding the formation of elliptical galaxies in clusters, the color-magnitude 
relation (OMR) plays a fundamental role. The CMR is an important tool, since in clusters it 
is tight and as a result provides a relation that can be used to tie together ages, metallicities 
and masses of galaxies. The scatter in the CMR is known to be small (~ 0.05 mag in the 
relation between V and U—V (see Bower, Lucey & EUis 1992a,b; Terlevich et al. 2001 for low 
redshift galaxies), and is also established in clusters at high redshift (Stanford, Eisenhardt 
& Dickinson 1998; Blakeslee et al. 2003), supporting the monohthic collapse scenario. The 
origin of the CMR has been explored extensively. A model based on the monohthic collapse 
scenario with a galactic wind has been able to explain the CMR as a sequence of metallicity 
(Arimoto & Yoshii 1987; Kodama & Arimoto 1997). This is not the only way to explain 
the CMR: models based on hierarchical clustering (Kauffmann & Chariot 1998; Ferreras, 
Chariot & Silk 1999; Terlevich et al. 1999; Ferreras & Silk 2000) producing a combination 
of age and metallicity variations could also explain the CMR. This controversial situation is 
almost due to the age-metaUicity degeneracy. In short, the origin of the CMR is still under 
debate. 

Spectroscopic properties, such as absorption line strengths, can partially break the age- 
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metallicity degeneracy. Worthey (1994) showed that the H/3 index, defined in Burstein et al. 
(1984), is sensitive to the luminosity- weighted age of galaxies and less to their metallicity. If 
a galaxy contains ionized gas, however, its emission makes the B.(3 absorption line weaker, 
and therefore an uncertain correction is necessary. To partially get around this one can use 
H7 and B.6 indices, such as H7A, H5a, H7F and H5f (defined in Jones & Worthey 1995 and 
Worthey & Ottaviani 1997), which are less affected by emission than H/5. These aforemen- 
tioned indices are some of the Lick/IDS indices (Burstein et al. 1984, Worthey et al. 1994), 
which are measuring line strengths of several absorption lines at low resolution. Although 
maybe less than colors, these indices are still affected by the age-metaUicity degeneracy. 

Partially to overcome the problem of low resolution Rose (1984) investigated spectra of 
elliptical galaxies using higher resolution indices (~2 A). His indices, such as RS/Fe I and 
H7/4325, are defined as ratios of neighboring absorption fines. Jones & Worthey (1995) 
defined a new index, H'Jhr that would be more sensitive to age. That index, however, 
is strongly affected by velocity dispersion (a) broadening of galaxy spectra, since the band- 
passes of feature and pseudocontinua arc very narrow, and therefore age determination using 
this index are difficult and has uncertainty. Later, Vazdekis & Arimoto (1999) introduced a 
new age indicator, H70-, which is very successful at breaking the age- metallicity degeneracy 
and usable for a wide range of a. While the presence of non-solar abundance ratio in most 
early-type galaxies (e.g., Gonzalez 1993; Kuntschner 2000; Poggianti et al. 2001, etc.) in- 
troduces systematic errors in the derived ages if Lick indices are used, errors from H7o- are 
less affected by non-solar abundance; i.e. errors mainly come from Poissonian noise. 

Vazdekis et al. (2001a) applied Rj^ to 6 elliptical galaxies in the Virgo cluster and 
demonstrated that it is a useful age-determinator for elliptical galaxies. The derived ages 
are uniformly old (>10Gyr), except for one galaxy (NGC 4478) which has bluer color, and 
the total metalficities and [Mg/Fe] ratios are increasing with luminosity, color and velocity 
dispersion. They concluded that the CMR is driven mainly by the mean stellar metallicity. 
Although this is an important result, it needs to be strengthened, since the number of galaxies 
in Vazdekis et al. (2001a) is small, and the luminosity spans only 3.3 magnitudes. 

In this paper we present new data of 8 elfiptical galaxies in the Virgo cluster observed 
with the Subaru Telescope together with the 6 galaxies of Vazdekis et al. (2001a). For the 
combined sample of 14 galaxies we present the analysis of the line strengths and determine 
their ages in various ways. We describe the details of our sample selection and observations 
in Section 2. In Section 3, we describe the line index measurements and stellar population 
analysis based on the simple stellar population (SSP) models of Vazdekis (1999). In Section 
4, we present the measurement of ages and metallicities of our sample, and discuss the 
advantage of using H70. index and compare our results with previous studies. In Section 
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5 the conclusions are given. In a future paper (Yamada et al. 2006) we discuss age and 
metallicity trends along the CMR and other scahng relations. 

2. Observations and Sample Selection 

In this paper we present data taken with the Subaru and the William Herschel Telescope. 
Vazdekis et al. (2001a) presented the first results of the WHT long slit observations. We 
re-reduced the WHT data to check the dependence of data reduction technique on the index 
measurements and added 8 new galaxies observed with the Subaru Telescope. 

2.1. Subaru Observation 

The Subaru observations were done on April 17-18, 2002 and June 26, 2003. The spec- 
troscopic observations were performed with the long slit mode of the FOCAS spectrograph 
with the MIT CCD and the R300 grism in the 2nd order (Kashikawa et al. 2002). This 
instrumental setup provides a dispersion of 0.67 Apixel"^ and a spectral range of AA ~ 
3800A-5800A. We used a slit width of 0".4, giving a resolution of 2.0 A (cr ~ 60 km s"^) on 
April 17, 2002 and a slit width of 0".6, corresponding to a resolution of 3.1 A (cr ~ 95 km s^^) 
on April 18, 2002 and June 26, 2003. The spatial scale was 0".21 pix^^ after binning 2 pixels 
and the seeing was 0".6-l".0 for all nights. Multiple exposures of 20 minutes length were 
acquired for each galaxy, as well as ThAr arc lamp exposures for wavelength calibration at 
every galaxy pointing. We obtained a number of dome and twilight flat fields for flatfielding. 
We also observed one spectrophotomeric standard star per run, Feige 34 in 2002 April and 
BD+33 2642 in 2003 June. 



2.2. WHT Observation 

The WHT observations were done in April 21-22, 1999. We used the ISIS double-beam 
spectrograph with the EEV12 CCD and R600B grating in the blue channel. The dispersion 
and spectral resolution is 0.44 Apixel"^ and 2.4 A (cr ~ 65 km s~^) provided by a sht width 
of 1".6 and a spectral range of AA ~ 4000-5500A. The seeing was 2"-3" during the first 
night and ~ 1".5 during the second. Multiple exporsures of 35 minutes length were obtained 
for each galaxy. We made use of tungsten fiat fields obtained at each galaxy pointing and 
twilight fiat fields as well as CuArNe arc lamp frames for A calibration for each exposure. 
For correction to a relative fiux scale the spectrophotometric standard star G191B29 was 
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used. 



2.3. Data Reduction 

Data reduction was done with the standard IRAF packages in the standard way; i.e. 
overscan correction, bias subtraction, flat flelding, A cahbration, sky subtraction and flux 
cahbration. We used dome flat flelds (for the Subaru data) and tungsten lamp flat fields 
(for the WHT data) for correction for pixel-to-pixel variations of the sensitivity. The spectra 
of WHT have wiggles coming from problems with the dichroic. We successfully removed 
these patterns by using tungsten lamp flat fields obtained for each galaxy pointing. We 
achieved AA (error in A) smaller than 0.07A for the A calibration, which is very important 
for the determination of H70-, since even a small shift in A can affect the K'-fcr measurement 
significantly (Vazdekis & Arimoto 1999). Cosmic rays were removed using the "cleanest" 
task of the REDUCEME package (Cardiel et al. 1998). Exposures whose spectra around H7 
feature were hit by cosmic ray were excluded. Flux calibration was done with the observed 
flux standard stars G191B29 (for the WHT data), Feige 34 and BD+33 2642 (for the Subaru 
data) . 

In this study, since we are interested in the central part of galaxies, we have summed 
up the spectra within rg/lO. Most of elliptical galaxies are known to have color and line 
strength gradients (e.g., Peletier et al. 1990, Gonzalez 1993; Kobayashi & Arimoto 1999). 
Here we used effective radii from the RC3 catalogue (de Vaucouleurs et al. 1991). Before 
summing the central galaxy spectra, we carefully aligned them with the rows of the array, 
and each row in spatial direction was corrected for the velocity shift due to the rotation by 
cross-correlating with the Vazdekis (1999) models. Variations of a within r^/lO were small 
(< 25 km s~^) for all galaxies, and therefore no correction for this effect was applied. 

To check for systematic effects due to telescope characteristics and data reduction, we 
compared the line indices for NGC 5831, a galaxy observed during both Subaru and WHT 
runs. We found excellent agreement between these indices measurements (see the open 
triangle (WHT) and open diamond (Subaru) in the middle panel of Fig. 2). Moreover, 
WHT data reduced independently by two authors (Y.Yamada. and A.Vazdekis) show also 
good agreement. 
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2.4. Sample Selection 

Table 1 lists the basic parameters of the observed galaxies, such as radial velocity (Vr), 
effective radius (re), velocity dispersion (a), visual magnitude {my), U — V and B — V colors, 
telescope, exposure time, slit width, and signal-to-noise ratio per Angstrom around 4340A. 

Our sample was primarily based on the Virgo galaxies observed by Michard (f982) and 
Bower et al. (f992a,b), which show the OMRs with surprisingly small scatter. We included 
3 elliptical galaxies from Jones (1998) in order to enlarge our luminosity coverage toward 
the fainter end. In Fig. 1 we reproduce the CMR of our sample galaxies. We selected the 
galaxies according to the following criteria: 1) Since we were only interested in elliptical 
galaxies, we excluded galaxies classified as SO in the references mentioned above and the 
RC3 catalogue (de Vaucouleurs et al. 1991). 2) We excluded galaxies having strong emission 
lines, such as NGC 4374 and NGC 4486. 3) We selected galaxies evenly for each luminosity 
bin niB < 11.5, 11.5 < < 12.5, 12.5 < rriB < 13.5, and itlb > 13.5, to cover a wide range 
of the himinosity. The luminosity spans 5.8 magnitudes, which is larger than the data of 
Vazdekis et al. (2001a), which spans only 3.3 magnitudes. 
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Fig. 1. — Color-magnitude relation of Virgo cluster elliptical galaxies (Bower et al. 1992a,b). Our sample 
is indicated by the large symbols with NGC numbers. NGC 4239, whose symbol is a filled triangle, has no 
U -V data. 
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2.5. Line Indices 

We used the Hj^- index (Vazdekis & Arimoto 1999; Vazdekis et al. 2001b) and the 
indices of the Lick system (Worthey 1994; Jones & Worthey 1995; Worthey & Ottaviani 
1997). For measuring the indices we used the FORTRAN program called "LECTOR" \ 
which is based on the equations (41)-(44) and simulations in Cardiel et al. (1998). 

We list the wavelength definition of 4 H'fcr indices in Table 2. Since the required S/N 
to distinguish between models of 12 Gyrs and 17.5 Gyrs is very high (S/N = 200 for 
H7o-<i30 and S/N = 400 for H7275), we integrated for a considerable time and conse- 
quently reached sound-to-noise ratios ranging from 115 to 490 for the selected apertures, 



See http: / / /www.iac.es/ galeria/ vazdekis/models.html 
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i.e., Te/lO. Although the S/N of our data is shghtly lower than hsted in Table 2 in most 
galaxies, we could achieve a high precision in the age determination for the following reasons: 
1) According to Cardiel et al. (2003), the real random errors of H70- are smaller than those 
computed by the LECTOR program since the wavelength definitions of H70- are not usual, 
i.e., the pseudo-continua and the feature bandpasswd overlap. We used the real error by 
using their new estimate of the random errors for H'-fa- 2) The 117^ index decreases rapidly 
for young ages, just as other Balmer lines do. Although the S/N A"-*^ for NGC4489 is low 
(~ 115), the error in age is only 0.6 Gyrs. 

The spectral range of our data goes from 3800A to 5800A to cover many of the Lick 
indices: B.6a and H7A (Jones & Worthey 1995), B.6f and H7F (Worthey & Ottaviani 1997), 
H/?, CNi, CN2, Ca4227, G4300, Fe4383, Ca4455, Fe4531, C4668, Fe5015, Mgi, Mg2, Mg6, 
Fe5270, Fe5335, Fe5406, Fe5709 and NaD (Worthey et al. 1994). Worthey (1994), Jones & 
Worthey (1995) and Worthey & Ottaviani (1997) have investigated the age and metaUicity 
sensitivity of these indices. They concluded that Balmer line indices are sensitive to age 
more than metaUicity. Besides H^y^ we used H/9 and B.Sf as age indicators for reference. 
Concerning the metaUicity, we used [MgFe] (defined in Gonzalez 1993) to estimate the mean 
metaUicity [M/H], 

[MgFe] = v^Mgfc • (Fe5270 + Fe5335)/2 . 

The feasibility of using [MgFe] for the total metaUicity is shown in Vazdckis et al. (2001a), 
Thomas, Maraston & Bender (2003), and Bruzual & Chariot (2003). Apart from this, we 
used several indices to determine abundances of several elements and molecules. For Fe we 
use the Fe3 index defined in Kuntschner (2000), 

Fe3 = (Fe4383 + Fe5270 + Fe5335) /3 . 

to give an Fe abundance estimate [Zpe/H]. Mg6, Ca4227, CN2 and G4300 are used to de- 
termine [ZMg/H], [Zca/H], [Zcn/H] and [Zch/H], respectively Tripicco & BeU (1995) and 
Trager et al. (1998) confirmed the dominant element or molecule in each line index. But each 
index is not pure, because the indices are still somewhat contaminated by other elements or 
molecules. For the G-band, Worthey (1994) and Tripicco & Bell (1995) show us that the 
contributcrs to the G-band are CH, Fe and Ti. Yet, since the main contributer is still CH 
(C), we write [Zch/H] for convenience. 

Note that, although the bandpasses are mostly those of the Lick system, we did not 
use the Lick/IDS system itself, instead we used fiux-calibrated spectra and no offsets were 
applied to reproduce the stars of the Lick/IDS system. Measured line strength of the indices 
are listed in Table 3 & 4. 
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0.030 


0.050 


0.065 


-2.202 


0.154 


0.202 


1.438 


5.649 


5.499 


0.013 


0.001 


0.001 


0.009 


0.016 


0.024 


-1.943 


0.116 


0.159 


1.701 


5.598 


5.672 


0.02:-! 


0.1)01 


D.OOi 


0,01!) 


0.():-!l 


0,011 



^Velocity dispersion measured from spectra. 

''Adopted total velocity dispersion for the measurements. 

<= Measurements for fTmeasured = 130 km s"! (NGC 4239, NGC 4387, NGC 4458, NGC 4467 and NGC 4489). 
^For H7(7 5 we show the measurement only for the appropriate velocity dispersion (Tmeasured 

(150, 190 or 300 km s"!). 



Table 4. Line Indices (Red) 



Galaxy Ca4455 Fe4531 C4668 H/3 Fe5015 Mgi Mg2 Mgb Fe5270 Fe5335 Fe5406 Fe5709 Fe5782 NaD 
Error (A) (A) (A) (A) (A) (mag) (mag) (A) (A) (A) (A) (A) (A) (A) 



NGC4239 


1, 


,612 


3, 


,456 


4, 


,839 


2, 


,204 


5, 


,497 


0, 


,031 


0, 


,150 


3, 


,135 


2, 


,856 


2, 


,641 


1, 


,746 


1, 


,029 


0, 


,819 


1, 


,793 


± 


0, 


,050 


0, 


,073 


0, 


,109 


0, 


,043 


0, 


,096 


0, 


,001 


0, 


,001 


0, 


,048 


0, 


,054 


0, 


,063 


0, 


,047 


0. 


,042 


0, 


,042 


0, 


,061 


NGC4339 


1, 


,818 


6, 


,097 


0, 


,098 


0, 


,251 


4, 


,495 


3, 


,200 


3, 


,065 


2, 


,090 


1, 


,064 


1, 


.015 


4, 


.229 


1, 


,785 


3, 


,584 


7, 


,859 


± 


0. 


,036 


0, 


,052 


0, 


,075 


0, 


,031 


0, 


,067 


0, 


,001 


0, 


,001 


0, 


,033 


0, 


.037 


0, 


.042 


0, 


.032 


0, 


,028 


0, 


,028 


0, 


,040 


NGC4365 


1, 


,544 


4, 


,090 


8, 


,207 


1, 


,467 


5, 


,670 


0, 


,137 


0, 


,301 


4, 


,810 


2, 


.886 


2, 


.552 


1, 


.754 




- 




- 




- 


± 


0, 


,026 


0, 


,039 


0, 


,060 


0, 


,025 


0, 


,057 


0, 


,001 


0, 


,001 


0, 


,027 


0, 


.031 


0, 


.035 





.026 




- 




- 




- 


NGC4387 


1, 


,703 


3, 


,632 


5, 


,927 


1, 


,662 


5, 


,897 


0, 


,096 


0, 


,234 


4, 


,200 


3, 


.141 


2, 


.964 


2, 


.004 




- 




- 




- 


± 


0, 


,049 


0, 


,074 


0, 


,116 


0, 


,048 


0, 


,105 


0, 


,001 


0, 


,001 


0, 


,051 


0, 


.057 


0, 


.065 


0, 


.049 




- 




- 




- 


NGC4458 


1, 


,678 


5, 


,080 


0, 


,068 


0, 


,206 


4, 


,241 


2, 


,799 


2, 


,607 


1, 


,671 


0, 


.966 


0, 


.897 


2, 


.686 


1, 


,513 


3, 


,392 


5, 


,040 


± 


0, 


,037 


0, 


,054 


0, 


,080 


0, 


,032 


0, 


,070 


0, 


,001 


0, 


,001 


0, 


,034 


0, 


.039 


0, 


.045 


0, 


.034 


0, 


,030 


0, 


,029 


0, 


,042 


NGC4464 


1, 


,609 


3, 


,458 


5, 


,103 


1, 


,599 


5, 


,332 


0, 


,104 


0, 


,248 


4, 


,425 


2, 


.829 


2, 


.506 


1, 


.745 




- 




- 




- 


± 


0. 


,039 


0. 


,059 


0, 


,093 


0, 


,038 


0. 


,084 


0, 


,001 


0, 


,001 


0, 


,040 


0, 


.046 


0, 


.053 


0, 


.039 




- 




- 




- 


NGC4467 


1, 


,733 


5, 


,620 


0, 


,099 


0, 


,243 


4, 


,234 


3, 


,121 


2, 


,897 


2, 


,022 


1, 


.026 





.918 


3 


.233 


1, 


,712 


3, 


,482 


5, 


,909 


± 


0, 


,042 


0, 


,062 


0, 


,092 


0, 


,037 


0, 


,082 


0, 


,001 


0, 


,001 


0, 


,040 





.045 





.052 





.039 


0, 


,034 


0, 


,034 


0, 


,048 


NGC4472 


1, 


,580 


5, 


,246 


0, 


,124 


0, 


,289 


4, 


,545 


2, 


,749 


2, 


,375 


1, 


,607 


0, 


.827 


0, 


.812 


4, 


.854 


1, 


,317 


3, 


,474 


8, 


,169 


± 


0, 


,036 


0, 


,052 


0, 


,074 


0, 


,030 


0, 


,067 


0, 


,001 


0, 


,001 


0, 


,034 


0, 


.037 


0, 


.042 


0, 


.032 


0, 


,029 


0, 


,028 


0, 


,040 


NGC4473 


1, 


,727 


3, 


,650 


8, 


,711 


1, 


,632 


5, 


,906 


0, 


,136 


0, 


,294 


4, 


,884 


3, 


.182 


2, 


.886 


1, 


.986 














± 


0, 


,022 


0, 


,034 


0, 


,052 


0, 


,022 


0, 


,049 


0, 


,001 


0, 


,001 


0, 


,024 


0, 


.026 


0, 


.030 





.023 














NGC4478 


1, 


,856 


3, 


,733 


6, 


,568 


1, 


,919 


6, 


,220 


0, 


,104 


0, 


,253 


4, 


,435 


3, 


.244 


3, 


.117 


2, 


.075 














± 


0, 


,045 


0, 


,068 


0, 


,106 


0, 


,043 


0, 


,095 


0, 


,001 


0, 


,001 


0, 


,046 


0, 


.051 


0, 


.058 


0, 


.043 














NGC4489 


2, 


,413 


6, 


,171 


0, 


,046 


0, 


,173 


3, 


,456 


3, 


,305 


3, 


,068 


1, 


,973 


1, 


.141 


1, 


.006 


2, 


.656 


1, 


,678 


3, 


,626 


6, 


,345 


± 


0, 


,051 


0, 


,074 


0, 


,108 


0, 


,042 


0, 


,093 


0, 


,001 


0, 


,001 


0, 


,047 


0, 


.052 


0, 


.060 


0, 


.046 


0, 


,041 


0, 


,040 


0, 


,059 


NGC4551 


1. 


,896 


6, 


,203 


0, 


,091 


0, 


,245 


4, 


,527 


3, 


,371 


3, 


,260 


2, 


,177 


1, 


.156 


1, 


.055 


3, 


.468 


1, 


,809 


3, 


,719 


7, 


,538 


± 


0, 


,033 


0, 


,048 


0, 


,071 


0, 


,029 


0, 


,063 


0, 


,001 


0, 


,001 


0, 


,031 


0, 


.035 


0, 


.040 


0, 


.030 


0, 


,026 


0, 


,026 


0, 


,038 


NGC4621 


1, 


,447 


3, 


,593 


8, 


,958 


1, 


,480 


5, 


,560 


0, 


,165 


0, 


,325 


4, 


,909 


2, 


.943 


2, 


.645 


1, 


.783 














± 


0, 


,013 


0, 


,019 


0, 


,029 


0, 


,012 


0, 


,028 


0, 


,000 


0, 


,000 


0, 


,013 


0, 


.015 


0, 


.017 





.013 














NGC4697 


1, 


,603 


5, 


,869 


0, 


,108 


0, 


,261 


4, 


,591 


3, 


,224 


3, 


,003 


2, 


,037 


1, 


.019 





.999 


4, 


.220 


1, 


,643 


3, 


,709 


8, 


,194 


± 


0, 


,020 


0, 


,030 


0, 


,043 


0, 


,018 


0, 


,039 


0, 


,000 


0, 


,000 


0, 


,019 


0, 


.021 


0, 


.025 


0, 


.019 


0, 


,016 


0, 


,016 


0, 


,023 



2.6. Stellar Population Models 



To derive the ages and metallicities of the galaxies, we used the revised evolutionary stel- 
lar population synthesis model of Vazdekis (1999), which has been made using the empirical 
stellar spectral library of Jones (1998) and updated with the new Padova isochrones (Girardi 
et al. 2000). The model predicts spectral energy distributions (SEDs) in the optical wave- 
length range, AA3976A-4456A and AA4764A-5460A, for simple stellar populations (SSPs) at 
a resolution of l.SA (FWHM). Previous models (e.g., Worthey 1994; Vazdekis et al. 1996) 
mostly used the Lick polynomial fitting functions for absorption line indices (Worthey et al. 
1994; Worthey & Ottaviani 1997) which are based on low resolution spectra (FWHM ~ 9A), 
thus limiting the application of the models to strong features. The new model provides flux- 
calibrated, high-resolution spectra, therefore allowing an analysis of galaxy spectra at the 
natural resolution given by their internal velocity broadening and instrumental resolution. 
Thus no uncertain correction for velocity broadening has to be apphed as is the case when 
using the Lick/IDS system. 



3. Age and Metallicity of Elliptical Galaxies 

Traditionally, a combination of metallic indices and Balmer line indices has been used to 
measure the ages of the central regions of observed galaxies (e.g. Gonzalez 1993; Worthey & 
Ottaviani 1997). Here we use a number of different age indicators: H70., H/3, RSp as defined 
by Vazdekis & Arimoto (1999), Worthey et al. (1994), and Worthey & Ottaviani (1997), 
respectively. For simplicity, we have divided our sample into three groups containing galaxy 
spectra with similar velocity dispersions and smoothed the model SEDs of Vazdekis (1999) 
to match the obtained total velocity dispersion. We chose (1) (Tmcasurcd ~ 150 km s"-*^, (2) 
(Treasured ~ 190 km s^^, and (3) (^measured ~ 300 km s^^. To allow direct comparisons between 
the galaxies having similar velocity dispersions, some small cr corrections were applied : the 
galaxies with (Xtotai (= (o"^aiaxv + '^instr)^) < ^-measured wcrc couvolvcd with the appropriate 
Gaussian to reach a fixed o"measured- On those spectra we determined our line strength indices. 
For example, the spectrum of NGC 4387 with a (Jtotai of 120 km s~^ was convolved with a 

Gaussian of (Jextra = 90 km S""^ to reach a simulated (Jmeasured (= ((^|alaxy + '^instr + ^extra)^) 

of 150 km s-^ 
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Fig. 2. — H7cr, H/?, HSp vs. [MgFc] diagrams. Top, middle and bottom panels are for galaxies with 
CTmoasurcd = 150 km s^^, 190 km s^^ and 300 km s^^, respectively. Model grids with various ages and 
metallicities are overplotted. Symbols are the same as in Fig. 1. 
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Fig. 4.- 



The Same as Fig. 3, but for H/3. Symbols are the same as in Fig. 1. 
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V 



3 3 4 

Fig. 5.— H7<,<i3o VS. [MgFc] diagram. Galaxies with a < 130 km s^^ (NGC 4239, NGC 4387, NGC 4458, 
NGC 4467 and NGC 4489) are plotted. The symbols are the same as in Fig. 1. 

3.1. H7(^ and Other Age Indicators 

In Fig. 2, H70. (left), H/3 (center) and RSp (right) are used as the age indicators, while 
the often-used metallicity indicator [MgFe] (Gonzalez 1993) is used as the metaUicity index. 
The velocity dispersion of the galaxies increases in each column from top to bottom. Note 
how the model grids change as a function of velocity dispersion. 

The left panel (H7CT-[MgFe]) of Fig. 2 confirms again that 117^ is well-suited to break 
the age-metallicity degeneracy (Vazdekis & Arimoto 1999). The lines of constant age are 
almost horizontal, and the grid of revised Vazdekis (1999) models is more orthogonal than 
the other two index-index diagrams (center and right panels) . Moreover, since B.'-/^ has much 
narrower index difinition than the other H7 indices such as H7A, metallicity and ^-enhanced 
effects should be significantly smaller (Thomas, Maraston & Korn 2004). Therefore, we are 
able to derive ages using H7(^ under small influence of metallicity, or abundance ratios such 
as [Mg/Fe], though it can be still affected by underlying isochrone structure with non-solar 
abundance ratios. Since H7 is less affected by emission lines (see Worthey & Ottaviani 1997, 
for example), the results from R'j^ are more reliable in general than those from H/3. 

Not shown here are diagrams using the other age indicators of Worthey & Ottaviani 
(1997), i.e., H7F, H7A, and H5a, as their abilities to disentangle ages and metallicities are 
similar or even smaller than RSp. In Fig. 3, we show 117^ as a function of several other metal 




-16- 



indices: Fe3, Mg6, Ca4227, CN2 and G4300. We also present the H/3 vs. metallicity indicator 
diagrams in Fig. 4 to demonstrate the disadvantage of using non-orthogonal diagrams for 
measuring ages of galaxies with non-solar abundance ratios. One can clearly see that the age 
obtained using H/3-Fe3 is older than that from H/3-Mg6, providing H/3 provide less accurate 
abundances of various elements. 

For galaxies with crtotai < 130 km/s (NGC 4239, NGC 4387, NGC 4467 and NGC 4489), 
the H7o-<i3o index can also be applied (Fig. 5). The H7o-<i3o index has a slightly different 
wavelength definition; i.e., a wider bandpass for the H7 feature than is the case for the other 
H7cr indices. It turns out that the ages derived from H7o-<i3o for these four galaxies are 
completely consistent with those from H7125. 

Balmer lines may be filled in with gas emission, becoming weaker and making it difficult 
to measure accurate mean ages. Although Gonzalez (1993) and Trager et al. (2000a) made 
corrections to the H/3 measurements using the [OIII]5007A emission line, we did not make 
any corrections here in this paper. We have investigated whether this correction is necessary 
by fitting the whole galaxy spectra with the SSP models to the data. If there is a residual at 
H/3 or at H7, it could be attributed to the presence of nebular emission. From these residual 
we found spectra that a only few galaxies show emission lines [OIII]4959 and/or [OIII]5007 
(NGC 4239, NGC 4489 and NGC 4697, sec Fig. 9). However they are small enough to be 
ignored for the H/3 and H7 age determinations. 

3.2. Ages from 117^ and H/3 

For each galaxy we determined the ages from R'jc^ and H/3 by comparing the data with 
the SSP models. The resulting ages are given in Table 5. The errors in the R'ja ages are 
dominated mainly by photon noise, while those for H/3 determined in this way are mostly 
due to systematic errors coming from the resulting non- orthogonal index-index model grids. 
Therefore, while the errors in ages from H70- can be improved by longer exposures, this is 
not straightforward for H/? ages. Note that in principle one can obtain better ages from H/? 
using more sophisticated modeling. In Table 5 we show the errors which one gets when using 
these simple index-index diagrams. We found that in most cases the youngest ages come 
from the H/3-Mg6 diagram, while the oldest ages come from H/3-Fe3 diagram. The errors in 
the H/3 ages fisted in Table 5 come from these two diagrams. One can clearly see that the 
observed galaxies cover a wide range of ages, from 3 Gyrs to over 15 Gyrs, particularly for 
the galaxies with lower velocity dispersion. 

Fig. 6 shows the age differences between H/3 and H'fa- They are consistent with each 
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other within the errors, since H/3 is not strongly affected by emission for our sample. Since 
[MgFe], used for deriving the B.p age, is almost insensitive to the [a/Fe] ratio, we infer that 
there is no systematic difference between the two ages, therefore, the B.f3 ages are also used 
in the following discussion. We note also that for several velocity dispersions many index 
- H/3 diagrams are still somewhat perpendicular, so that also the systematic errors are not 
much larger than the errors from photon noise. 



Fig. 6. — Comparison between ages from H70. and H/3. Symbols are the same as in Fig. 1. 

Fig. 7 gives the comparison between our results and the literature, i.e., Gorgas et al. 
(1997), Trager et al. (2000a), Terlevich & Forbes (2002), Proctor & Sansom (2002) and 
Caldwell, Rose & Concannon (2003). Gorgas et al. (1997), Trager et al. (2000a) and Proctor 
& Sansom (2002) measured ages from H/3 using Worthey (1994) 's models. In Trager et al. 
(2000a), H/3 measurements were corrected for emission using the [OIII] 5007A hue. Caldwell 
et al. (2003) used a combination of Balmer lines (Hn/Fe) and their new SSP models. The ages 
of Terlevich & Forbes (2002) arc a compilation of several line index papers (Davies et al. 1993; 
Carollo et al. 1993; Fisher ct al. 1996; Gonzalez 1993; Goudfrooij et al. 1999; Halhday 1998; 
Kuntschncr 1998; Longhctti ct al. 1998; Mchlcrt ct al. 1997; Vazdekis & Arimoto 1999). All 
age measurements except for those of Vazdekis & Arimoto (1999) come from B./3. Our relative 
ages are consistent with those of the previous papers, although our results are systematically 
older. The systematic difference can be explained in several ways: 1) over correction in 
the literature for emission line contamination on H/3. To correct the emission effect in H/3 
measurements, all the past studies used the AH/3 = 0.7[OIII] equation of Gonzalez (1993), 




AflB (Gjr) 
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which had been derived from the residuals between galaxy spectra and stellar spectra (not 
SSP models). We made no correction for emission lines, since our measurements are based 
on the H7 feature, which is almost free of emission lines due to residual star formation or 
planetary nebulae. 2) Worthey (1994) 's model, which the previous works used, is based 
on polynomial fitting functions of the Lick/IDS system. To compare our results with the 
previous literature, it is necessary to convert our measurements of the indices to the Lick/IDS 
system. Having done that, we find no systematic difference between our measurements and 
the literature measurements which use the Lick/IDS system, although the error is large 
due to our relatively small number of reference stars used for the conversion. Accordingly, 
we conclude that the systematic difference in ages does not come from differences in the 
observational configurations used nor from data reduction errors. The models of Vazdckis 
(1999) employs the Padova Isochrones (Girardi ct al. 2000), while the Worthey (1994) 
model used the Revised Yale Isochrones (Green et al. 1987). Since the Padova Isochrones 
are systematically hotter, this may explain the difference in the estimated ages. 3) The age 
estimates corresponding to the H/5 vs. [MgFe] diagram in the Lick/IDS system differ from 
those estimated using flux-calibrated response curves at the speciflc velocity dispersion of 
the galaxy (see Kuntschner et al. 2002). 

We note that NGC 4365 has a very old age, of more than 20 Gyr (but consistent with 
14 Gyr within la). We need S/N/A> 350 for galaxies with a > 250 km s~^, while we have 
achieved in this case S/N/A= 250. As a result the uncertainty in the age is rather large. 
Note that the H/3 age of NGC 4365 is 14.0 Gyrs, which is consistent with that of Davies et 
al. (2001). Finally note that relative age differences are much more accurate than absolute 
ages, due to possible calibration problems in the stellar evolutionary tracks (Vazdekis et al. 
2001b). 
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Fig. 7. — Comparison between our ages and ages from previous studies. T2000, TF2002, C2003, G1997 
and PS2002 indicate Trager et al. (2000a), Terlevich & Forbes (2002), Caldwell et al. (2003), Gorgas et al. 
(1997) and Proctor & Sansom (2002), respectively. Symbols are the same as in Fig. 1. 
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Table 5: Age and Metallicity of Virgo Elliptical Galaxies 



NGC Age(Gyr)" 


[M/H]" 


Age(Gyr)'' 


[M/H]'' 


IZFe/H]<= 




lZca/H]'= 


[ZcN/ar [Zch/h]'' 












[ZMg/ZFel° 


[Zca/ZFel° 


[ZcN/ZFel° [ZcH/Zpel° 



4239 


6 2-1* 


„ ,-.-0.11 
-0.16+0 JO 


6 0-1° 
"■"+0.5 


n 1C-0.05 
-"■10+0.05 


-0 14-°°'* 

^■^^+0.09 


-o^is^S.i" 
-o^o4;o°.SJ 


+o.2o;°;^i 

+o^34;S:?5 


„ ,-.-0.07 
-0.16+0.08 

-0 02-° °l 
o-O'^+o.oo 


-0 32-° l° 
"■•3-^+0.16 

„ ,0-0.10 
-"■i»+0.07 


4339 




-0 01-° °^ 
"■"^+0.08 


8 9-° '' 


+0 09-° °* 

-t-u.uy+0.04 


-0 12-°°'' 
"■^■'+0.06 


+0 17-°13 
+"-i'+0.09 

+0 29-°-°'' 


-0 21-°11 
"■''i+o.ii 

-0 09-°-°^ 
"■""+0.05 


_l_n 1K-O.O6 
+0.16+0 06 

, „ oR-O.OO 
+ U.^» + 00 


-0.68+0 16 

e.c-0.06 
-0^56+o.io 


4366 




+ U.IU+0 j4 




_i_n iK-° °<' 
+016+0.05 


-0 13-°°'^ 

"■^■^+0.07 


+0.35+0 20 

+o^48;S,-J^ 


,.,.-0.09 
-"■»f'+0.16 

07-0-02 
+ "■"'+0.09 


+ 20-°-ll'' 
+ "■"^" + 0,09 

-0.00 
+ "■•'•' + 0.02 


-"■('i;S.o5 
+o^i2;S:oo 


4387 




-0 12-"°^ 
"I'^+O.ll 


i2-9;!:I 


-0 12-°°^ 
°l^+0.07 


-0 14-°°° 

"■^*+0.08 


-o^o2;S:!2 
+0 12-°°^ 

+U. 1^+0 04 


-0 20-° °^ 
"■^^"+0.13 

-0 06-° °° 
"■""+0.05 


-0 03-" °'^ 

"■"^+0.07 

+ 11-001 
+ "■11+0.04 


-0 14-°^3* 
"■1^+0.21 

+o-oo;S:?3 


4468 


13 7-^ '' 


-0 24-° °® 

"■^^+0.08 


i2.8;J:^ 


-0 20-° °^ 
U-^U+0.03 


-0 37-0 °3 
°-^'+0.07 


-0 03-° °° 

"■"•^+0.13 

+0 34-° °3 
-t-u.oi+o 06 


n oK-0.06 
U.OD+0 06 

+0 01-° °^ 
-t-u.ui+o 00 


+o^o8;S:o°J 
+o^45;S:So 


-0.69+0.21 

-0 32-°-°'' 
" -^'=+0.14 


4464 


18-5+Ji% 


-0 30-° °'' 
u.ou+0.08 




-0 21-° °=' 
"■''^+0.03 


„ qo-0.03 


-0 04-° °® 
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" Metallicities and pseudo-abundance ratios come from the ti-y^ vs. Fe3, Mgb, Ca4227, CN2 and G4300 diagram, respectively. 
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3.3. Metallicity and Abundance Ratios 

By breaking the age- metallicity degeneracy, we can estimate the mean metallicity [M/H] 
and a heuristic abundance ratio [Zx/Fe] for the elements or molecules X that happen to 
dominate the absorption of the indices plotted in Fig. 3. 

[M/H] is determined from the H7o--[MgFe] and the H/3-[MgFe] diagram. Thomas et al. 
(2003) have shown that the [MgFe] index is a good indicator of the total metallicity [M/H] 
and almost independent of [a/Fe] (see also Vazdekis et al. 2001b; Bruzual & Chariot 2003). 

Abundance ratios contain important information about the formation of galaxies, and 
are understood to indicate different time-scales for various types of Supernovae (SNe). The 
enhancement of a elements by SNe II explosions occurs without delay from the onset of star 
formation because their progenitors are massive stars. On the other hand, SNe la, whose 
progenitors are white dwarfs of binaries, mainly supply iron peak elements 1-1.5 Gyrs after 
the onset of star formation. As a result, lower [Mg/Fe] indicates that the contribution of SNe 
la is larger; i.e. star formation is relatively longer (Worthey 1998) or the number fraction of 
SNe la progenitors is enhanced. Concering CN, both carbon and nitrogen are ejected from 
AGB stars whose life time are placed between those of SNe II and SNe la. Therefore CN 
could be another potential indicator for the duration of star formation (Sanchez-Blazquez 
et al. 2003), although the photospheric abundance of carbon and nitrogen are rather fragile 
to the convective dredge-up in the stellar interior and thus may not directly reflect the 
nucleosynthesis history. However, it has been recently found a correlation between the 
abundance ratio of CN with respect to Mg and Fe and X-ray luminosity of galaxy cluster 
(Carretero et al. 2004). 

Since the Worthey (1994) models only have solar abundance ratios, past studies (Trager 
et al. 2000b; Proctor & Sansom 2002) used the response functions on the Lick/IDS system 
of Tripicco & Bell (1995) to derive abundance ratio of each element. Since the method 
using their response functions is complex and introduces unknown systematic errors, we 
did not use it. We therefore adopted a more straightforward method using the H7o--metal 
index diagrams, i.e., we directly estimated [Zpc/H], [Zng/H], [Zca/H], [Zcn/H] and [Zch/H] 
from the H7cr vs. Fe3, Mg6, Ca4227, CN2 and G4300 diagram with the SSP model grids, 
respectively. For galaxies outside the model grid, we slightly extrapolated the grid to measure 
their metallicities. 

Although our method is more straightforward, the estimated abundances still include 
uncertainties. For example, it is well known that a line index has contamination from other 
elements. Thomas et al. (2003) have pointed out the [a/Fe] bias of stellar libraries using 
stars in our Galaxies, i.e., metal poor stars have positive [a/Fe]. Another problem is that 
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the SSP models used here have been made using stellar interior models and isochrones with 
solar abundance ratios. One therefore has to keep this in mind, although it is safe to use 
them for deriving relative ages as we have shown in this work. To differentiate between our 
heuristically estimated abundances and true [X/H] abundances, estimated abundances in 
this paper are written as [Zx/H] and so on. 

In Fig. 3 we confirm the well-known result that elliptical galaxies have non-solar abun- 
dance ratios (e.g. Worthey et al. 1992; Trager et al. 2000a). For example, NGC 4365, 
NGC 4472, NGC 4473, NGC 4551, NGC 4621 and NGC 4697 show a large enhancement in 
[Zms/H] and many of our sample galaxies, except for NGC4239, NGC4387 and NGC4489 
show significant enhancement in [Zcn/H]. [Zca/H] is subsolar for most galaxies, which was 
first pointed out by Vazdekis et al. (1997), for reasons that are not fully clear yet. One would 
expect Ca to behave like other ct-elements such as Mg, on the basis of a nucleosynthetic and 
chemical evolution sense. Vazdekis et al. (2003) showed that Call triplet indices in the 
infrared are sensitive the global metallicity for [Fe/H] < —0.5 and depend very much on the 
IMF for [Fe/H] > -0.5. 



3.4. Fitting with SSP Models 

We show the observed and the best matching SSP spectra in Figs. 8 & 9 ^. Both the 
observed spectra and the SSPs are divided by a fitted polynomial of order 2 and normalized 
at 4350A and 4950A for the blue and red spectra, respectively. We do not fit wavelength 
shorter than 4000A since no data are available there (WHT) or they have low sensitivity 
(Subaru). The observed galaxies are indicated by thick lines, and models by thin lines. Ages, 
metallicities and o"totai of the models are quoted in each figure. We also present the ratio of 
observed spectrum to the SSP model in the lower panel of each diagram. 

For NGC 4489, we examined the models with several ages and metallicities to confirm 
the accuracy of our age and metallicity estimate from the line indices (Fig. 8). Since the 
estimated age and metallicity of this galaxy arc 3.1 Gyrs and [M/H] = +0.20, respectively, 
we show the fits with age of 2.00, 3.16, 5.01 Gyrs, and with [M/H] of +0.20, ±0.0, -0.38. 
For reference, we also show the cross correlation peak heights derived using the IRAF task 
"fxcor", whose value of 1.00 means that the perfect fit has been reached. Fig. 8 shows that 
the model of the same age and metallicity derived by H'fa, i.e., 3.1 Gyrs and [M/H]=+0.20, 
clearly gives the best fit. 



^Fig. 9 are available in the electronic version. 
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Fig. 8. — The full spectral fitting for NGC 4489 using the SSP model spectra of various ages and metallic- 
ities. Thick and thin lines indicate spectra of galaxies and SSP models, respectively. The ratio between the 
observed spectra and the best fit model is shown in the lower panels. Note that we keep the same scale for 
residuals. 
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Fig. 9. — Spectral fitting for galaxies in the upper panels. The observed spectra and the best fit model 
spectra are indicated by thick lines and thin lines, respectively. We present the ratio between the observed 
spectra and the model spectra in the lower panels. We note [M/H], Age and cxtotai of the model in the figure. 
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We have confirmed that most of galaxies show excellent fits to models of ages and 
metallicities derived from the H7o--[MgFe] diagram. For NGC 4365, whose Hj^- ^-ge is over 
20 Gyrs, we estimated that the age is 14-18 Gyrs by fitting of models with various ages. One 
can see the small emission lines of [0111]4959A and/or [Olll]5007A for NGC 4239, NGC 4489 
and NGC 4697, but we cannot see residuals in either H7 or H/? for any of the galaxies. Hence, 
we conclude that emission lines do not affect our age measurements. The overabundance of 
CN can also be seen in the fits. The observed spectra of NGC 4365, NGC 4458, NGC 4464, 
NGC 4467, NGC 4472, NGC 4621 and NGC 4697 are below the models around the CN 
molecular band at ~ 4150A. NGC 4464, NGC 4472 and NGC 4621 show stronger absorption 
of the Mg triplet line at ~5180A. In conclusion, our age and metaUicity measurements agree 
with the whole spectral fitting and they can be trusted in our discussion. 

4. Conclusions 

We have obtained very high S/N (>100, per A) spectra of 8 elliptical galaxies selected 
along the CMR in the Virgo cluster with the Subaru Telescope and 6 eUiptical galaxies 
observed with the WHT previously reported in Vazdekis et al. (2001a), spanning a total 
luminosity range of 5.8 magnitude. By using the H7(^ method, we have determined ages, 
metallicities, and abundance ratios of 14 elliptical galaxies in the Virgo cluster. 

The R'ja index is almost independent of the galaxy metaUicity; i.e., it can break the 
age-metallicity degeneracy, therefore the age from the H70. index does not depend much on 
metaUicity ([M/H]) or abundance ratios. Although, like all indices, H'y^ is not completely 
immune from abundance ratios (Vazdekis et al. 2001b), the relative ages remain secure. 
The resulting ages show a large spread, from 3 Gyrs to over 15 Gyrs. We have compared 
these ages with previous studies that used the H/3 index and find that the relative ages are 
consistent. 

We find that Zca (from the Ca4227 index) is following Zpc, and that Mg (from the 
Mgb index) is overabundant with respect to Fe. These results are partially known in past 
studies, however, previous approaches were affected by the age-metallicity degeneracy. If 
one takes those abundances at face value, it is not easy to reconcile underabundant Ca with 
overabundant Mg at the same time, since both Ca and Mg are thought to come from SNell. 
Our results seem to suggest that Ca and Mg are not produced in SNell of exactly the same 
mass range. 

We also fitted the spectra with SSP models finding good agreement with our age and 
metaUicity determination. Contamination of emission in H/3 and H7 is small enough not to 
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affect our age measurements, while three galaxies (NGC 4239, NGC 4489 and NGC 4697) 
have emission line in [OIII]4959 and/or [OIII]5007. Our fits show the overabundance of some 
elements in several galaxies: NGC 4365, NGC 4458, NGC 4464, NGC 4467, NGC 4472, 
NGC 4621 and NGC 4697 have a strong CN band at ~4150A and NGC 4464, NGC 4472 
and NGC 4621 show stronger absorption in the Mg triplet than the fitted SSP model. 

In a subsequent paper, based on the ages and metallicities determined here and together 
with the velocity dispersions and luminosities, we will discuss the star formation history of 
elliptical galaxies along the CMR of the Virgo cluster. 
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